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c^ ■ Abstract 

We describe a simple lattice model of higher-curvature quantum gravity in two 
dimensions and study the phase structure of the theory as a function of the 
curvature coupling. It is shown that the ensemble of flat graphs is entropically 
unstable to the formation of baby universes. In these simplified models the 
growth in graphs exhibits a branched polymer behaviour in the phase directly 
before the flattening transition. 



1. Introduction 

A longstanding problem of quantum gravity is to determine the effects of the addition 
of higher-curvature counterterms to the canonical gravity action. This problem is most 
tractable in two- dimensions, mainly because the theory can be regularized and studied 
using the discretization approach (see [jl| for reviews). In this case, the simplest term 
of this type that one could add is KR{gY, where R{g) is the scalar curvature of a two- 
dimensional metric g and the coupling constant n is inversely proportional to the square 
of the worldsheet ultraviolet cutoff. The continuum Euclidean partition function is given 
by the path integral over all worldsheet metrics 



Z = jDgexp - j^cfz ^d^g (^^ + ^-^R{g) + KR{g)' 



:i-i) 



where /i is the cosmo logical constant and G is the gravitational constant. For k ^ the 
R^ term in the action is irrelevant and the system lies in the pure gravitational phase, i.e. 
that with action consisting of the cosmological and (topological) Einstein terms and with 
critical string exponent 7str = —\ (for a fixed spherical topology). On the other hand, 
the limit k ^ oo suppresses large curvature fluctuations of the metric and the statistical 
ensemble becomes flatter and flatter at short distance scales. At /t = oo only flat surfaces 
contribute to the partition function. 

The problem of whether or not this system undergoes a phase transition between the 
flat and pure gravitational phases at some finite coupling k^. has been studied using quan- 
tum Liouville theory in and investigated numerically in [Q. An exact non-perturbative 
solution has been obtained in |Q for spherical topologies using matrix model techniques 
and it was shown that 7str = —\ for all finite values of k. The model always reduces 
at large length scales to a model of pure gravity, i.e. there is no transition to a non- 
perturbative phase of flat metrics, and only when all surfaces with non-minimal R{gY 
are completely forbidden (k = oo) does 7str change. The key feature of this proof is an 
appropriate extended model that interpolates between the fixed lattice (or crystalline) 
flat space model and the pure gravity model p. A simplified model which exhibits the 
same characteristic features has been studied in 0. In this paper we shall give a sim- 
ple demonstration of the absence of a phase transition from a random to a flat phase of 
R^ gravity directly in terms of the fractal structure of two-dimensional quantum gravity. 
This has the advantage of being technically much simpler than the solution given in 
while at the same time exposing some of the physical and geometrical characteristics of 
R? gravity. 

In section 2 we shall describe how to incorporate the effects of the R^ action in two- 
dimensions, using the fact that the value of 7str is controlled by the structure of baby 
universe formation on the surfaces of the ensemble. This leads to the study of surfaces with 
branching outgrowths (which we call "bottleneck surfaces" ) which have a contribution to 
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the R^ term proportional to the neck thickness. In section 3 we prove that the generating 
function for the bottleneck surfaces is an analytic function of the curvature coupling 
K, thereby demonstrating the absence of a phase transition as k, is continuously varied. 
Finally, in section 4 we present a model which captures the essential qualitative features 
of the flattening transition. There we prove that the ensemble of fiat graphs is unstable to 
the formation of baby universe outgrowths, so that the flattening transition can only take 
place exactly at k = oo. These results show explicitly that no matter how fiat the system 
is at short distance scales, it always destabilizes at long wavelengths into the familiar 
ensemble of highly fractal baby universes. 

2. Baby Universes and Lattice R^ Coupling 

The dynamically triangulated version of pure two-dimensional quantum gravity is 
given by the partition function 

m = Ee-^^^(^) = E^''^ E 1 (2-1) 

A A TsTa 

where Z{A) is the fixed area partition function and Ta is the topological class of trian- 
gulations of area A {A is proportional to the number of triangles of 7^). A given choice 
of T G 7^ corresponds to a particular discretization of the surface. The local intrinsic 
curvature of a vertex i G T of coordination number g^ is Ri = 27r(6 — qi)/qi, so that the 
discretized form of the R^ action is given by 



{qi - 6) 
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d'z Jdet g R{gf -^ An' E ^^^-^^ (2-2) 
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However, there is a much simpler way to incorporate the effects of the R' term by 
simply finding an appropriate extended model which interpolates between pure two- 
dimensional quantum gravity and a fiat phase of worldsheet metrics. The main idea 
comes from the fact that the universal constant 7str, which describes the nature of the 
geometry, is related to a surface roughness structure of two-dimensional quantum gravity, 
the distribution of "minimal bottleneck baby universes" on the ensemble of triangulations 
[[7|. A baby universe of area Aq is a small region of the triangulation joined to the bulk 
solely by a minimal neck, i.e. a loop consisting of only three links. Intuitively, they can 
be pictured as forming bubblings on the surface. For a triangulation of area A — ;> oo (the 
continuum limit), the asymptotic behaviour of the fixed area partition function is 

Z{A) ~ A^^'^-^ e'^"^ (2.3) 

where /xq is the critical cosmological constant (so that the series ( ^.1[ ) converges for all 
/i > /io)- From this it can be shown [0] that the distribution r2^(ylo) of baby universes is 
related to the string susceptibihty exponent by 

nA(Ao) ~ kA^-^^'^{A - Aoy^'^-'Al^'^-' (2.4) 



where here and in the following k denotes a constant which is independent of the areas 
and is of the order of unity. This result shows that, as 7str increases, the number of 
baby universes increases as nA{Ao)/A ~ A;Ao"'"~ . Thus 7str not only measures the fractal 
structure of the typical surfaces that contribute to the partition function, but also their 
branching ratio into minimal bottleneck baby universes. 

Now consider the situation where the curvature coupling k is very large and the / R^ 
term is thereby trying to make the surfaces as smooth as possible. A surface of spherical 
topology has a minimum J R^, denoted -Rj^j^, which is independent of its area (basically 
corresponding to a dual graph consisting of hexagons and 12 pentagons). If we take two 
such surfaces (not necessarily of the same area) and join them as shown in fig. 1 by cutting 
out a triangle from each one and stitching them together, then we create a surface with 
one minimum length bottleneck baby universe and J R^ = 2/?^;^ + A, where A is the 
contribution to the / R^ term arising from the bottleneck. We can iterate this procedure 
to produce surfaces with n bottlenecks and 



R' 



-'^min ~r " l-'^min 



A 



(2.5) 



The important point here is that / R^ is not proportional to the area, but only to the 
number of bottlenecks. On the other hand, most of the surfaces of a given area created in 
this way are geometrically very different from the single surface with minimum J R^. We 
can extend this exercise by joining surfaces together with larger, non-minimal necks of 
length i. In this case it is clear that the contribution to / R^ can be at most proportional 
to the length i of the bottleneck. 




(u) 



(r) 
^„(A-Aq) 

Figure 1: Gluing an unmarked baby graph of area Aq (counted by VIq (Aq)) onto a rooted 
n-bottleneck graph of area A — Aq (counted by Vt'^^A — Aq)) along a common deleted disk 
(shaded areas) to produce a rooted (n + 1) -bottleneck graph of area A — B — C . 



To get an idea of how the inclusion of the R^ term modifies the distribution functions, 
it is instructive to consider a simple example. Consider the ensemble of closed rooted trees 



on a fixed lattice. The total number of such trees of length l (taken to be the number of 
links in the tree) is J\f^^\£) which behaves asymptotically as (see for example Pj) 

ArW(£)~r° e"^ (2.6) 

where Oq is a universal constant which is independent of the particular type of lattice, 
while p > is a non- universal constant. For a given tree of length i there are generically 
£ different ways of drawing the tree on the lattice. We introduce an extra fugacity factor 
Q-Ki £qj, each loop, which is the analog of the curvature term J R^. The total number of 
closed rooted trees on a fixed lattice of area A is then 

W{A,k) = J2iAf^'\i) e""' (2.7) 

e 

For large areas A -^ oo, we can replace the sum in ( p.7|) by an integral over i G [0, oo) 
and suppress the trees which grow too large {i"^ > A) by an exponential damping factor 
g-£ M. Thus asymptotically we have 

W(A, k) - r di (iU^''\f.) e-"^ e-^'/^ ~ itA"«+3/2 ^{p-nfA/i (2.8) 

We see that the only effect of the fugacity coupling is to renormalize the entropy parameter 
as p — *> p—n. The distribution function (p.8|) is analytic in k and the inclusion of a fugacity 
factor for the loops is irrelevant from this point of view. In particular, we may adjust the 
coupling so that k = p, in which case VV(v4, p) has a power law growth with the area. Thus 
the non-universal tree growth has no effect on the critical behaviour of the corresponding 
grand canonical ensemble quantities either. Note that for a self-avoiding random walk on 
a regular square lattice, we have Oq = — | (see for example 0), and W(A, k) contains a 
purely exponential growth with the area A of the lattice. By tuning to the point n = p, 
the distribution in this case is independent of the size of the lattice. 



3. Analytic Structure of the Bottleneck Graph Generating 
Function 

We will now write down an expression for the partition function of bottleneck surfaces 
and argue that it is an analytic function of the R^ coupling n. We work at fixed spherical 
topology and consider the generating function for the ensemble of rooted bottleneck graphs 
which is given by 

fiM(p,«:;^,W}) =Ee-'^^E^"fiM(A«:;W}) =Ee-'^^E^" E ^b (3.1) 

An An BeB^ 

where B'^ denotes the set of all n-bottleneck surfaces of area A. The local weight of a 
graph B G ;B^ is given by 

^B = n «<?.-! e""'" (3.2) 

6eB 



where b are the bottleneck nodes of B of branching number g^ and neck thickness ib, 
and am are arbitrary weights associated with the nodes. The number Q^'^ (A) of rooted 
spherical graphs of area A with n bottlenecks (where the area of each triangle is taken 
to be unity) can be constructed inductively as follows. For n = this number is related 
to the number VIq (A) = Z{A) of unmarked spherical graphs as VIq (A) = AQ.Q (A), 
since there are A possible triangles to mark in an unmarked graph. To construct a rooted 
graph with n + 1 bottlenecks, we cut out a disk of area B and perimeter i from a rooted 
graph with n bottlenecks, and a disk of area C with the same boundary length i from a 
0-neck unmarked graph (fig. 1). Gluing the two cut graphs together along the common 
perimeter of the deleted disks yields a rooted graph with n + 1 bottlenecks. In doing so, 
one must be careful of the degeneracies which occur in this inductive cutting procedure. 
Consider two spherical (0-neck) graphs which are identical to each other except for their 
triangulations inside two disks, of areas Bq and Cq and with identical boundary length i, 
on each respective graph. These two inequivalent graphs, when glued onto other lattices 
to produce bottleneck graphs, yield the same set of surfaces. 

We shall first consider the sub-ensemble of linear chain bottleneck graphs. We intro- 
duce the disk amplitude Zi{i; B) which is defined to be the number of surfaces of area B 
and one boundary of length i. If we join two surfaces of areas Bq and B — Bq > Bq, each 
having a single boundary of length i, along their common boundaries, then we obtain a 
closed surface of area B with a marked loop of length i partitioning it into two parts of 
areas Bq and B — Bq. Since the boundary has i links, there are generically i different 
ways of joining the two surfaces to obtain distinct final surfaces (provided that Bq and 
B — Bq are large enough). Moreover, any surface of area B with a marked loop of length 
C. partitioning it into two parts of areas Bq and B — Bq can be uniquely constructed in 
this way. Therefore, for large B and Bq., there are 

Gi{B, Bq- i) ~ e Z,{£; BQ)Z,{i; B - Bq) (3.3) 

closed surfaces of area B with a marked loop of length i that partitions the surface into 
two parts of areas Bq and B — Bq. The asymptotic large area behaviour of the disk 
amplitude is 

Zi{i; B) ~ BT='^-2 e^^"-^ r-f^tr-i ^pi (34) 



where p > is as in ( |2.6| ). This amplitude therefore increases with i for small i, and for 
large perimeter loops (£^ > B), Zi{i] B) is exponentially damped by worldsheet finite size 
cutoff factors like e-^'/^ 0. 

We shall also use the annulus amplitude Z2{£i,£2', B) which is the number of surfaces of 
area B with two holes of boundary lengths Ci and C,2- Its large area asymptotic behaviour 
is 

^2(^1, ^2; B) ^ B^^'^-^ e^«^ (^1 + ^2)"'^=*^"^ e''(^i+^2) (3.5) 

and for large loops it is suppressed by the worldsheet infrared cutoff q^^^^+^^) I^ g. Then 



there are 

G2{B, Bi, B2; £1, £2) = £x£2 Zi{£v, B,)Z,{£2; B2)Z2{£i, £2; B - B, - B2) (3.6) 

distinct surfaces of area B with two non-intersecting loops of lengths £1 and £2 enclosing 
areas Bi and -B2, respectively. Neglecting the graphs counted by Zi{£; B) and Z2{£i, £2', B) 
which themselves have bottlenecks, i.e. assuming that (|3.4|) and (|3.5| ) correctly count the 
disk and annulus graphs with no bottlenecks of lengths smaller than £ and £1, £2, an ra-neck 
linear chain can then be constructed by gluing together n — 1 annulus graphs in between 
two disk graphs. We assume that there is a one-to-one correspondence between the latter 
configurations of graphs and those obtained by slicing an n-bottleneck graph along each 
of its necks. These assumptions will be sufficient to determine the analytic properties of 
the full generating function. 

We may now write down an expression for the number C^^\A;£i^ . . . ,£„) of rooted 
linear graphs of total area A and with n bottlenecks of lengths £1 , . . . , £„. In the continuum 
limit, we can replace sums over areas and bottleneck loop lengths by integrals. The lower 
limits of integration for the perimeter integrals are 3 links (the minimal bottleneck size), 
while those of the area integrals are 1 triangle. As all area integrations are ultraviolet 
finite, there is no danger in continuing their limits down to Ai = 0. Incorporating the 
smooth infrared cutoffs on the lengths as described above and the loop fugacity factor 
g-2K£i for each neck, this yields the combinatorial identity 



/ra+l „oo \ 

£M(A,«:;£i,...,C) = ny dAA 6{A-J:^Ai)A,£^Z^{£,■A,)e-^^'^e-'y^^ 



xZi(4;A„+i)e-^"/^"+^ 



i=\ 



It is convenient to analyse this relation in the grand canonical ensemble. For this, we 
define the Laplace transforms of the loop amplitudes by 



Zi(£;/i) = / dB e-^"^ Zi{£; B) e- 
Jo 



e/B 



(yZ/i - /JQ ) 
£2 



l-7str 



eP'K^^,^,\2^/jr^^^. 



Z2(£i,£2;/i) = / d5e-^^Z2(£i,£2;5)e" 
Jo 



{h+i2Y/B 



(V/^-/^o)' 



-7str 



,P(4+^.) K_^^^^ [2./jr^o {£1 + £2)) (3.8) 



where we have used ( p.4|) and (|3.5| ), and Ky is the irregular modified Bessel function of 
order v. The corresponding expression in the grand canonical ensemble is thus 

^ /"OO 

'-'n [f^j l^'-j £l} ■ ■ ■ }£n) = / ciA e L,^ [A, k; li, . . . ,ln) 

Jo 
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(V/^ - /^o 



,l-(n+l)7str 



'l^n 



n-1 



X 



e2(p-«)E.^« -Q 



i=l "-« 



^-7=tr (2V/^-/^0^0^1-7str (2V^-^0C 



ir_,,,,(2v/;?r^7^(£, + £,+i)) (3.9) 



+ 



-j+i 



We are now interested in the analytic properties of the total distribution of linear 
n-bottleneck graphs 






(3.10) 



as a function of k G R"*". Because of the ultraviolet cutoff on the length integrations, 
possible singularities in (|3.10|) would come only from the behaviour of the integral as 
£j — i> oo. For our purposes we may therefore approximate the modified Bessel functions 
in (|3.9|) by their asymptotic behaviours K,^{z) ~ z~^^'^ e~^ for |z| — ;> cxo and consider the 
function 



^^:\f^,^) 



(v/;?i^7^)-^^"-'^-^"+'^"='^ 



(v/2) 



n+l 



X 



n 

Vi=l 



n-1 






(3.11) 



In order that (|3.11|) be generically convergent, we must be in the phase with fi > fi, 
where 

{p - Kf 



/ic = /LtO + 



(3.12) 



At n = He there is a phase transition, but the critical behaviour is just that of the usual 
continuum limit of the discretized surface model. This will follow from the fact that both 
(|3.11 ) and ( p.l2|) are analytic functions of k, and the usual critical point pc = Po can be 
reached by tuning the R^ coupling to the value k, = p. We now focus on the analytic 
properties of C^\p, k) assuming that p > pc- It can be uniformly bounded by using 
the inequality ti + fj+i > ^J2liti^ to write down an upper bound on the integrations on 



the right-hand side of (|3.11| ) which can be evaluated in terms of the incomplete gamma- 
function r(a;,x) = /^(itt"~^e~* to give a bound on the derivatives of C^\p, k) with 
respect to k, 

i(n-l)-(n+l)7st 



d'ci:Hp,K) 



dn^ 



< 



( V/^ - /^o )' 



2 a//! — /io + K — p 



-(n-5)/2+fen 



2(7n-ll)/4-fc(n+l) 

xr(-i +k,Q{2^jr^Q + K~p))T{-\ +k,Q{2^jr^Q + K-p) 

(3.13) 



r ( i + A;, Q{2^p - po + K- p 



n-2 



for /c > 0. 

For non-integer or negative a, r(a,x) is a multi- valued function of x with a branch 
cut along the negative real x-axis. In the phase with p > p^, the linear n-bottleneck 



distribution is therefore uniformly bounded by an analytic function of 1^ [i — /io + i^ — p- 
In terms of the generating function for the linear bottleneck graph ensemble, defined as 
a formal power series 

oo 

£«(/i,«:;;2) = ^z"£W(/i,«:) (3.14) 



n=0 



in a variable z, this bound reads 



d''B'-\^i,K-. 



dn^ 



< 



r(-i +fc,6(2v/^r^7^ + fi:-p))r(-| +k,6{2./]r^o + K-p 



■)2fc+3 



r Q +k, Q{2y/p - po + K- p) 
r Q + /c, 6(2v//i - Aio + /€ - p) 



-l/2+k 



(3.15) 



Being uniformly bounded by an analytic function, the generating function cannot have 
any singularities as k > p—2^p — po is varied and it is analytic in the variable 2^ p — /io + 
K — p. Furthermore, since the critical cosmological constant ( |3.12| ) is itself an analytic 
function of k, there are no singularities as k is varied throughout its range. Therefore, 
the linear bottleneck generating function is an analytic function of k G M^. There are 
no phase transitions within the ensemble of linear bottleneck graphs as one continuously 
varies the B? coupling constant. 

We now consider the full bottleneck graph ensemble with generating function ( |3.1| ). 
To study a bottleneck distribution with arbitrary branchings, we shall need to use the 
genus n-loop amplitude ZniJ-x-, . . . ,in', B). We can determine its approximate asymptotic 
behaviour as follows. For each additional loop that is drawn on a graph there is an extra 
combinatorial area factor giving an overall entropy B"Z{B) for locating the centers of 
the loops on the closed surface. Furthermore, each loop on the surface contributes an 
exponential length growth e''^% and the smooth infrared cutoff on large perimeter loops 
is e~^^i^'' '^. The crucial point, however, is that the n-loop amplitude depends only on 
the sum of the loop lengths ^ . The natural ansatz is then 

/ n \ -7str-l 



Zn{i 



1; • • • ) ■tni 



B] 



^-fstr-3+n gMoB 1^2^ 



^PJ2^ 



-(E,^-)'/^ 



(3.16) 



u=l 



We note that this ansatz only dictates the large area and loop length dependence of the 
amplitude, but not its coefficient which depends on n. 

To establish (|3.16|) , we use induction on n and the consistency condition 

rB 

(^1 



^n-ll^l; 



,L^i;B) = J^ dBo J^ d£Z^{i,,...,L_,,i;Bo)Z^{t,B-Bo) (3.17) 

which represents the combinatorics of gluing together disk amplitudes and ra-loop ampli- 
tudes to generate (n — l)-loop graphs. Using the ansatz (p.l6|) and changing variables in 
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the area integral to x = Bq/{B — Bq), the right-hand side of (|3.17|) becomes 



where L„_i = X^iL/^i- For S ^ oo we may approximate the integral over x in ( 3.18|) by 
its asymptotic behaviour which integrates to the modified Bessel function 



+ Ln-l 



7str — 1 



K,^^^^,{2i{i + L^^^)/B) 



(3.19) 



From the asymptotic behaviour Kv{x) ^ T^ ^T{i')x '^ + . . . for x — > and z/ 7^ 0, it 
follows that the large-area behaviour of ( p.l8 ) is thus 






V-l^t.+l ^£ ^ ^^_^^_T,,,,_i ^2pl-2l{i+L„^,)/B ^3^20) 



Within the present approximations the loop integral (|3.20|) may be evaluated using the 
saddle-point approximation. The stationary condition is 



2(p 



L 



n-l 



£(£ + L 



4f{i + Ln-i) 



n-l 



B 



(7str + 1)^ - (7str - 1)(^ + ^n-l) = (3.21) 



Solving the cubic equation (|3.21|) and evaluating the one-loop fiuctuation integral cor- 
responding to (|3.20| ) for the root which is positive and regular at Ln^i = 0, we find in 
the limit B -^ oo that ( ^.20| ) is proportional to the left-hand side of ( ^.17| ) as given by 



(|3.16|) . Taking the Laplace transform of the expression (|3.16| ) , we find the corresponding 
amplitude in the grand canonical ensemble, 



Zniil, ■ ■ ■ Jn] fJ') - ( V/^ - /^O ) 
~ (V/W - /iO 



2-n-7str 



2-n-7st 



■ft 



n— 3 



g(p-2v/^:=7:^)^^4 for I 



oo 



(3.22) 



With the bottleneck node couplings a^ completely arbitrary, the full generating func- 
tion may be generated by the recursive equation which is represented symbolically in fig. 
2. Summing over all configurations shown there gives the identity 

dAi dA2 6{A-Ai-A2) rdiie-^^z['~\i;Ai)Q(^i,A2,K;z,{am}) 



oo /n+1 „oo \ / " /-c 

AZ{A) + Y.aJl[ dA] 6{A-E^aMI[ 

n 



dL 



Kti 



(3.23) 



where Z^\ii, . . . ,in; A) = A Zn{ii, . . . ,in', A) is the marked n-loop amphtude and we 
have introduced the generating function Q{i, A, n; z, {am}) for bottleneck surfaces of area 
A and with a single deleted disk of boundary length i. The left-hand side of ( p.23|) is 
the desired quantity fl^^''{A, k; z, {am}) and from ( p.22| ) it follows that the corresponding 
relationship in the grand canonical ensemble is 



i1('")(/i, k;2;, {«„} 



" di i e^P-^^' K_^^^^ (2^jr^^i) n[£, /i, k; z, {am}) 

OO 1 _. _ 

(/i - flo)^"^"'" + J2'^n (V/^ - /^O ) 
n=l 

n l°° de, £, e^P-""^'^ n[e,, /i, k; z, {am}) 




n-2 



Kl-n-u, (2 V/i - /iO Ei (-i) 



(3.24) 



This iterative relation defines a complicated integral equation for VL{(i, /i, /?; z, {0;^}) whose 
solution gives the generating function (|3.1|) . The complexity of the identity ( p.24| ) makes 
an explicit solution or even an analyticity analysis intractable. However, one may argue 
that the solution to ( 3.24|) is an analytic function of k as follows. A generic branching 
distribution function for the bottleneck ensemble will always involve products of functions 
of the form ( |3.22|) . From the inequalities 




< 



n 



(n-l)/n 



n. 



„(n-|)/n 



n < A 
n> A 



(3.25) 



it follows that a general bottleneck function can always be bounded from above by a 
function which is given by a product of incomplete gamma-functions and other elementary 
functions. In the phase /x > /Xc, the n-neck distribution functions are analytic in p. The 
grand canonical distribution function f2*^^)(/i, k; z) is therefore expected to be an analytic 
function of k G M"*". Note that the recursive definition (13.241) is reminescent of that for 



the generating function of a branched polymer ensemble |jTO[. In the next section we 
shall consider a slight simplification of the model defined by ( |3.24|) which is amenable to 
explicit analysis and thereby demonstrate that this similarity is not a coincidence. 



4. A Model for the Flattening Transition 

We shall now present a model for the transition between fiat and random graphs which 
captures the essential qualitative features of the flattening transition. We will consider 
two simplifications of the full bottleneck graph generating function. First, we start with 
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Figure 2: Schematic representation of the iterative definition of the bottleneck surface 
generating function. The shaded circles represent the loop functions Q{£, A, k; z, {am}) 
while the lines represent 0-neck surfaces. Each attachment is done by gluing the surfaces 
together along a common deleted disk. 



the ensemble of regular square lattices,[| which in the schematic representation of fig. 2 are 
represented by the rooted lines. Let Z^^ (A) be the number of such flat graphs of area A 
and with a single square loop of side i drawn on them (the 0-bottleneck partition function) . 
We are interested in the generating function Ze{A) for the ensemble of bottleneck surfaces 
of area A with a single loop of side i drawn on them, which is depicted by the shaded 
circles in fig. 2. It is constructed recursively as described at the end of the previous 
section. We start with a flat base lattice of area Aq and put an outgrowth of bottleneck 
graphs with a square loop of side ra on it. There are Aq possible ways to do this and the 
area of the new graph is Aq — tl^. However, the effective area remaining to insert a second 
outgrowth of bottleneck surfaces of loop perimeter An is Aq — 4n^ because of the excluded 
area effect. Each such outgrowth is attached to the base graph with a fugacity e~'^'"("\ 
where w{n) is a function which satisfies w{l) = 1 and w'{n) > 0, but which does not grow 
as fast as n'^. 

This procedure can be used to build the full generating function Zi{A) by summing 
over all possible ways of attaching outgrowths to the base graph. Dividing by the permu- 
tation symmetry factor for a /c-branch surface, we have 



Ze{A,K) 



'(0) 



^V^ Ao/4n2 



zriA)+ E zriAo) E E i. 



Ao=r^ 



X 



n=l k=l 

oo 



k\ 



Ai=n2 



i=l 



6(^A-Ao + kn^ - Ei A 

(4.1) 



*Actually, one must consider square lattices with positive curvature defect insertions to be 
able to close them on a spherical topology. This modification would not affect the qualitative 
behaviour described in the following. The model could of course be extended to include non- 
square lattices and arbitrarily shaped bottlenecks. 
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It is again convenient to analyse the grand canonical partition function which here is 
defined by the discrete Laplace transform 



no o V^ 



A=^2 „=1 



A/in^ 



(4.2) 



where we have used the binomial theorem to do the sum over k from ( [4.1D . Upon carefully 
interchanging the sums over n and A in (|4.2|) we arrive at 



i,(/i, n) = nf^ (/x) + Y: Zf {Un) + E Z^^ i'^n) (4.3) 

n=l n=£ 

where we have introduced the function 

7^/) (/i) = if) (^) - - E v^ e-"^ ^f (^) ^ (4-4) 

and the rescaled cosmological constants 

/x„(/x, «:) = /i - ^ log (l + 4n2 e^"'-'^-^-) Z„(/x, «)) . (4.5) 

Note that the grand canonical base graph generating function behaves asymptotically for 
large loops as 

#)(/i)=e-^^^^,(/i) (4.6) 

where geifJ') < 5'i(/^) and gi{^) converges for all positive fi and diverges at // = 0. It 
follows that the sums in ([4.3|) converge provided /i„(/i, n) > 0. 

We will now analyse the analytic properties of the solution of the infinite-term recursion 
relation (|4.3|) as a function of /x and k. From ([4. 21) it follows that 



i,(/i,/t = oo) = if)(/x) (4.7) 

so that the suppression of curvature fiuctuations via an infinite i?^-coupling leaves only 
fiat (square lattice) graphs. In this fiat phase, the partition function diverges only at the 
critical point /i^°) = and formally the string susceptibility exponent is 7str = 2. For 
ft < oo, it follows immediately from ( [4.3|) that the partition function Zi{fi, k) determines 
a different universality class than the fiat phase. To see this, we suppose that Zi{fi, k) has 
the same critical behaviour as Z\ (f^), i.e. that Zi{fj., k) diverges at some critical point 
/i = ^c- If this were the case then the smallest of the /in(/i, i^) in (|4.5| ) would become zero 
at some fi = fi'^> fie, and so from ( |4.3|) and ([4.6|) the partition function Zi{fi, k) diverges 
before its critical point is reached. Thus a critical behaviour at any k, < oo whereby 
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Zi{ii, k) itself diverges is inconsistent. It is only at k = oo that the universality class of 
the random geometry changes. 

Physically then, the flat lattice system at k = cxd is entropically unstable to the 
formation of baby universes that make up the full gravitational ensemble, which gives a 
worldsheet geometric picture of why there is no flattening phase transition in the two- 
dimensional R^ quantum gravity model. In the remainder of this paper we will deduce 
what the nature is of the phase when one perturbs the fiat graph ensemble by baby 
universe outgrowths as described above. To get an idea of what the system at k < oo 
represents, consider the simplification whereby we allow outgrowths of only a single loop 
side i. This means that we keep only the n = C. term in ([4.1|) and proceeding as in the 
general case we would then arrive at 

Z,{^i,K)=zf\^^,). (4.8) 

This is the standard behaviour of the generating function for a branched polymer ensemble 
[pT| . Indeed, as before, Z^^fj,, k) cannot diverge at the critical point fic because then 
/2^(/i, k) would reach before ^^(/i, k) reaches cxd. Differentiating both sides of ( |4.8| ) with 
respect to /x using (|4.5|) yields 

?(0)/ 






(4.9) 



Since dZ^ {ftij/dpe -^ — cxo as j2e -^ 0, it follows that the denominator of ( [4. 91) van- 
ishes at some finite value of fie where the numerator is regular and non- vanishing. Thus 
dZilfi, K)/dfi diverges at the critical point, and generically we get a branched polymer 
ensemble with string susceptibility exponent 7str = +| |n[. However, the actual bottle- 



neck ensemble contains a complicated mixing of all loop lengths and more care must be 
exercised in deducing the critical behaviour. 

We begin by deducing some analytic properties of Z^^fi, k) for k < cxd using the 
recursion relation (|4.3|) . Note that for fixed i, Ze{fi, k) is a monotonic decreasing function 
of /i G M"*", while for fixed /x it is a decreasing function of C.. Let us examine the behaviour 
of the functions /i„ very close to the critical point /x = jic- Since /Xc = at k = oo, at very 
large k we expect that /ic is very small so that for ra^yU <^ 1 we find 



/^n+l ~ /^n > 



^ log (l + 4n2 e^"'-'^-(") Z„(/i, k)) 

-^^^2 log (l + 4(^ + 1)' e^("+^)'---("+^) i„(/i, n)) 



^log(l + 4n2e-'^'"(")Z„(^, 



4n2 

- ^^^\lY log (l + 4(n + If e-^'"^-+''> i„(/ie, k) 
> (4.10) 
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for fi r^ fic and k ^ 1. Thus {fin} is an increasing sequence for small n. Denoting 
/imin = infn/in, it follows froui (^4.2|) and ( [4.6|) that 



2i{ft,K) < e 



/^min-*- 



fi-ll/^r 



+ 



TT 



2 V A^n 



(4.11) 



On the other hand, because /imin < yU, from ( |4.3| ) it follows that for large i we have 

Z^(/i, k) > e"''"""^' (4.12) 

The bounds ( [4.11|) and ( |4.12| ) show that for large i we may write 



Zi{jJ,K) 



f^rain*-- 



/(^,/i) 



(4.13) 



where the function f{i,fi) is larger than 1 and grows at most linearly with i. We then 
see that at large enough n 



log4n^ log4(n + l)^ Kw{n + 1) 

H'n+l H'n — !l n 77 '. TTn ~r 



An? 



4(n + l) 



4(n + l) 



Kw[n) 
Ain? 



(4.14) 



so that now, depending upon the nature of the function w(n), the /i„ may start to fall 
again. However, we also have that 



3 1 

/ioo = lim /in = 7/W + T/Wmin > yWii 
n^oo 4 4 



(4.15) 



so that we can identify /imin = /^i- It follows from these properties of the sequence {fin} 
that the infinite series in ( |4.3| ) is absolutely convergent. 

We now rewrite (14.31) in the form 



(4.16) 



Z,(/x, k) = Zi(/x, /.) + ??/)(/.) - ??i°)(/.) + Y: {^Pifin) - Z^ifir 



n=l 



t7(o). 



The right-hand side of ( [4. 161 ) depends only on the known functions Ti.^ (/i) and the par- 



tition functions Zn{fi, k) ioi 1 < n < i. Iterating (|4.16|) thus determines Zn{fi, k,) as a 



function only of /i, k, and Zi, for any n > 1. Using ( ^4.3| ) we see that Zi satisfies 

Ziifi, k) = nf\fx) + T (/i, /£, Zi(/i, K)) (4.17) 

where implicitly the function T is given by 

oo 

T{fi,K,Z^) = E^^H/^n) (4.18) 

n=l 

finif^, K,Zn) = /i-^log(l + 4n2e^"'-'^'"(")i„) (4.19) 



together with ([4.16| ), but explicitly it is a very complicated function. Moreover, ( [4.19|) is 
an analytic function of fi, k, and Zi. 
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Differentiating (|4.171 ) with respect to /i gives 

dZA 



dfi 




1- 





(4.20) 



The numerator of ( [4.2(]| ) is a well behaved function for yU > but at some critical point 
the denominator will vanish and thereby generate a non-analyticity in Zi. To see this, 
we note from ([4.19|) that 

^] =y^-'«/„ (4.21) 

where we define 

^^'" = -ZP'ifiri) e"'(4'^"-3'^)-'^'"(") (4.22) 

and 

fe^^] = 1 + E (-4''" - ^"'i fn (4.23) 

C^l / ^l,K n=l 

where we have used ( [4.16|) . Note that A^'"' is a monotonically decreasing function of i, so 



that dg) implies 



/,+! - /, = E (^'''''" - ^''") /n < (4.24) 

n=l 

provided that /« > Vn > 1. Thus {/^} is a monotonically decreasing sequence. Under 
these same assumptions we also have 



n=l n=l 

^''"'^^^ E /n(^V(/in) - ^;(/in)) e^^'^'^'^-^-^-''^'^"). (4.25) 



< e 

ra=l 



The last line of ( [4. 25] ) vanishes in the limit ^ — » oo, so that the first series in the recursion 



relation ( [4.23| ) vanishes as i —>■ oo. Using this fact and (^4.23|) it follows that ( [4. 211) can 
be written in terms of the single quantity /oo = lim^^oo fi as 

7^ = 1-/00. (4.26) 

oZi 

Again assuming that the sequence {fi} is monotonically decreasing we have that the limit 
of the sequence has the lower bound 

00 
/oo>l-E^'''" (4-27) 

n=l 

which is positive at large enough /i. On the other hand, as /i decreases the sequence will 
go negative and then start to oscillate. It is clear that at small enough n this must happen 
because it is straightforward to choose a /x for which /2 is already negative. It follows 
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by continuity that the critical point where the sequence {fi} monotonically decreases to 
/oo = must exist, and therefore the system is a branched polymer. 

To determine the value of 7str we consider two cases. If w{n) is a strictly increasing 
function of n, then by working at very large k we can ensure that the first term in ( [4.21|) 
dominates so that at the critical point 



[91^1) - 9[mi)} e-'' (4.28) 

and higher terms are suppressed by factors of e"''''"''-"^""'^^''-'. As long as {fi} is a decreasing 



sequence, the sum of terms with n > 2 in ( [4.2 1|) converges to something negligible com- 



pared to the first term. Now we can make the standard argument that the denominator 
of (|4.2CI| ) will vanish linearly with Zi{fi, k) which we assume to behave as 

Zi(/i, k) = Zi(^/ic(/t),/«) -5(/ic(fi:),fi;) (/i-/ic(«:)) '*' (4.29) 

where the function B is regular at the critical point, and therefore 7str = |. In the case 
that w{n) is constant this argument fails. To show that the model has generic branched 
polymer behaviour we would need to show that 

^ when /oo = (4.30) 



dZl 



IM,K 



but we have not succeeded in finding a proof of this. However, it is straightforward to 
solve the recursion equations ( [4.23| ) numerically and fig. ^ shows the behaviour of /oo for 
e"'' = 10~^. It is clear that 7str = 2 ^^ ^^^^ '^^^^ ^^ well. 

We conclude that the system defines a generic branched polymer phase of two dimen- 
sional quantum gravity with string susceptibility exponent | as soon as the coupling k is 
tuned away from infinity. This model only allows a rather simple subset of all possible 
branchings and we expect that the ordinary pure gravity phase (with 7str = — |) is recov- 
ered by including all of the more complicated branchings in which the path on the base 
graph is allowed to be arbitrary. The basic point is that the model above gives a physical 
understanding for the absence of a fiat phase when the curvature coupling constant k, is 
large. 



This work was supported in part by PPARC grant # GR/L56565 and the Danish 
Natural Science Research Council. 
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log( M-u ) 



Figure 3: Logarithmic plot of f^o against fi for the case w{n) = 1, e " = 10 '^. The 
crosses represent computed values with fj,c = 0.040112 and the straight line has slope |. 
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